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c^ . Abstract 

In this paper, the hnear free flexural vibration of cracked functionally graded 

material plates is studied using the extended finite element method. A 4- 

p. • noded quadrilateral plate bending element based on field and edge consis- 

^ . tency requirement with 20 degrees of freedom per element is used for this 

Q^ ! study. The natural frequencies and mode shapes of simply supported and 

rn I clamped square and rectangular plates are computed as a function of gradi- 

C^ ■ ent index, crack length, crack orientation and crack location. The effect of 

thickness and influence of multiple cracks is also studied. 



Keywords: Mindlin plate theory, vibration, partition of unity methods, 
extended finite element method. 



1. Introduction 

Engineered materials such as laminated composites are widely used in auto- 
motive and aerospace industry due to their excellent strength-to and stiffness- 
to- weight ratios and their possibility of tailoring their properties in optimizing 
their structural response. But due to sudden change in material properties 
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between the layers in laminated composites, these materials suffer from pre- 
mature failure or by the decay of stiffness characteristics because of delam- 
inations and chemically unstable matrix and lamina adhesives. The emer- 
gence of functionally graded materials (FGMs) [H |2] has revolutionized the 
aerospace and aerocraft industry. The FGMs used initially as thermal bar- 
rier materials for aerospace structural applications and fusion reactors are 
now developed for general use as structural components in high temperature 
environments. FGMs are manufactured by combining metals and ceramics. 
These materials are inhomogeneous, in the sense that the material properties 
vary smoothly and continuously in one or more directions. FGMs combine 
the best properties of metals and ceramics and are strongly considered as a 
potential structural material candidates for certain class of aerospace struc- 
tures exposed to a high temperature environment. 

It is seen from the literature that the amount of work carried out on the 
vibration characteristics of FGMs is considerable [3llll[5l[6l|71IHll9]. B.e et 
al., [5] presented finite element formulation based on thin plate theory for 
the vibration control of FGM plate with integrated piezoelectric sensors and 
actuators under mechanical load whereas Liew et al., [Bjhave analyzed the 
active vibration control of plate subjected to a thermal gradient using shear 
deformation theory. Ng et al., [7] have investigated the parametric resonance 
of plates based on Hamilton's principle and the assumed mode technique. 
Yang and Shen [TD] have analyzed the dynamic response of thin FGM plates 
subjected to impulsive loads using Galerkin procedure coupled with modal 
superposition method whereas, by neglecting the heat conduction effect, such 
plates and panels in thermal environments have been examined based on 
shear deformation with temperature dependent material properties in [8]. 
Qian et al., [IT] studied the static deformation and vibration of FGM plates 
based on higher-order shear deformation theory using meshless local Petrov- 
Galerkin method. Matsunaga [4J presented analytical solutions for simply 
supported rectangular FGM plates based on second-order shear deformation 
plates. Vel and Batra [3] proposed three-dimensional solutions for vibrations 
of simply supported rectangular plates. Reddy [12] presented a finite element 
solution for the dynamic analysis of a FGM plate and Ferreira et al., [H] 
peformed dynamic analysis of FGM plate based on higher order shear and 
normal deformable plate theory using the meshless local Petrov-Galerkin 
method. Akbari et al., [13] studied two-dimensional wave propagation in 
functionally graded solids using the meshless local Petrov-Galerkin method. 
The above list is no way comprehensive and interested readers are referred 



to the literature. 

FGM plates or in general plate structures, may develop flaws during man- 
ufacturing or after they have been subjected to large cyclic loading. Hence 
it is important to understand the dynamic response of a FGM plate with 
an internal flaw. It is known that cracks or local defects affect the dynamic 
response of a structural member. This is because, the presence of the crack 
introduces local flexibility and anisotropy. Moreover the crack will open and 
close depending on the vibration amplitude. The vibration of cracked plates 
was studied as early as 1969 by Lynn and Kumbasar ^14j who used a Green's 
function approach. Later, in 1972, Stahl and Keer [15] studied the vibration 
of cracked rectangular plates using elasticity methods. The other numeri- 
cal methods that are used to study the dynamic response and instabihty of 
plates with cracks or local defects are: (1) Finite fourier series transform [T6] : 
(2) Rayleigh-Ritz Method [T7]; (3) harmonic balance method y^; and (4) 
finite element method [191 [20]. Recently, [21] proposed solutions for the vir- 
bations of side-cracked FGM thick plates based on Reddy third-order shear 
deformation theory using Ritz technique. Kitipornchai et al., [22] studied 
nonlinear vibration of edge cracked functionally graded Timoshenko beams 
using Ritz method. Yang et al., [23] studied the nonlinear dynamic response 
of a functionally graded plate with a through-width crack based on Reddy's 
third-order shear deformation theory using a Galerkin method. 
In this paper, we apply the extended finite element method (XFEM) to study 
the free fiexural vibrations of cracked FGM plates based on first order shear 
deformation theory. We carry out a parametric study on the infiuence of 
gradient index, crack location, crack length, crack orientation and thickness 
on the natural frequencies of FGM plates using the 4-noded shear fiexible 
element based on field and edge consistency approach [21] . The effect of 
boundary conditions and multiple cracks is also studied. Earlier, the XFEM 
has been applied to study the vibration of cracked isotropic plates [25 | l26 t |27] . 
Their study focussed on center and edge cracks with simply supported and 
clamped boundary conditions. 

The paper is organized as follows, the next section will give an introduction 
to FGM and a brief overview of Reissner-Mindlin plate theory. Section [3] 
illustrates the basic idea of XFEM as applicable to plates. Section Hj presents 
results for the free flexural vibration of cracked plates under different bound- 
ary conditions and aspect ratios, followed by concluding remarks in the last 
section. 



2. Theoretical Formulation 

2.1. Functionally Graded Material 

A functionally graded material (FGM) rectangular plate (length a, width 
b and thickness h), made by mixing two distinct material phases: a metal 
and ceramic is considered with coordinates x, y along the in-plane directions 
and z along the thickness direction (see Figure ([2])). The material on the 
top surface {z = h/2) of the plate is ceramic and is graded to metal at 
the bottom surface of the plate {z = —h/2) by a power law distribution. 
The homogenized material properties are computed using the Mori-Tanaka 
Scheme EHl 



Estimation of mechanical and thermal properties 

Based on the Mori-Tanaka homogenization method, the effective bulk mod- 
ulus K and shear modulus G of the FGM are evaluated as [2S1 123 EHl ITT] 
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Here, Vi (i = c, m) is the volume fraction of the phase material. The sub- 
scripts c and m refer to the ceramic and metal phases, respectively. The 
volume fractions of the ceramic and metal phases are related by Vc + Vm = 1, 
and Vc is expressed as 

K(-)=(^)", n>0 (3) 

where n in Equation (J3]) is the volume fraction exponent, also referred to as 
the gradient index. Figure ([1]) shows the variation of the volume fractions 
of ceramic and metal, respectively, in the thickness direction z for the FGM 
plate. The top surface is ceramic rich and the bottom surface is metal rich. 
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Figure 1: Through thickness variation of volume fraction 




Figure 2: Co-ordinate system of a rectangular FGM plate. 

The effective Young's modulus E and Poisson's ratio u can be computed 
from tlie following expressions: 
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The effective mass density p is given by the rule of mixtures as [3] 

P = PcVc + PmVm 



(4) 



(5) 



The material properties P that are temperature dependent can be written 
as BT] 



P = PoiP-iT-' + 1 + PiT + P2T' + P^T' 



(6) 



where Pq-, P-i, Pi, P2, P3 are the coefficients of temperature T and are unique 
to each constituent material phase. 

2.2. Plate formulation 

Using the Mindlin formulation, the displacements u, v,w at a point (x, y, z) in 
the plate (see Figure ([2])) from the medium surface are expressed as functions 
of the mid-plane displacements Uo, Vo, Wo and independent rotations 9x, Oy of 
the normal in yz and xz planes, respectively, as 



u{x,y,z,t) 
v{x,y,z,t) 
w{x,y,z,t) 



Uo{x,y,t) + z9^{x,y,t) 
Voix,y,t) + zey{x,y,t) 
Wo{x,y,i) 



(7) 



where t is the time. The strains in terms of mid-plane deformation can be 
written as 
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The midplane strains Sp, bending strain £f,, shear strain Sg in Equation 
are written as 
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where the subscript 'comma' represents the partial derivative with respect to 
the spatial coordinate succeeding it. The membrane stress resultants N and 
the bending stress resultants M can be related to the membrane strains, Sp 
and bending strains Sb through the following constitutive relations 
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where the matrices A = Ajj,B = Bij and Db = Dij] {i,j = 1,2,6) are the 
extensional, bending-extensional coupling and bending stiffness coefficients 
and are defined as 



rh/2 _ 

{A,, 5.,, A,}= / Q.,{1, Z, Z^] dz (11) 

J -h/2 

Similarly, the transverse shear force Q = {Qxz,Qyz} is related to the trans- 
verse shear strains Eg through the following equation 

Q,j = E.jSs (12) 

where Eij = /_;j /2 Q'^i'^j dz; {i,j = 4,5) is the transverse shear stiffness 
coefficient, Vi,Vj is the transverse shear coefficient for non- uniform shear 
strain distribution through the plate thickness. The stiffness coefficients Qj,- 
are defined as 
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where the modulus of elasticity E{z, T) and Poisson's ratio u are given by 
Equation (jlj). The strain energy function U is given by 

U{S) = ]- f {ejAsp + ejBe, + eJ^Sp + e^De,. + ejEs,} dfi (14) 

where 5 = {u.v^w^Ox^Oy} is the vector of the degree of freedom associated 
to the displacement field in a finite element discretization. Following the 
procedure given in [32], the strain energy function U given in Equation ^\M 
can be rewritten as 

U[5) = 1<5^K5 (15) 

where K is the linear stiffness matrix. The kinetic energy of the plate is 
given by 

ns) = \l [uiil + vl + wl) + h{el + el)] dn (le) 

t/ ill 



where /q = /_^/2P(-2) dz, Ii = X-L2 -^^P(-^) dz and p{z) is the mass density 
that varies through the thickness of the plate given by Equation ([5]). Substi- 
tuting Equation ( !T5|) - (TT6|) in Lagrange's equations of motion, the following 
governing equation is obtained 

M5 + K5 = (17) 

where M is the consistent mass matrix. After substituting the characteris- 
tic of the time function [33] 8 = —u'^S, the following algebraic equation is 
obtained 

(K - w^M) 6 = (18) 

where K is the stiffness matrix, u is the natural frequency. 

3. Overview of the extended finite element method 

In this section, we give a brief overview of the XFEM for plates. By exploit- 
ing the idea of partition of unity noted by Babuska et al, [M], Belytschko's 
group [35] introduced the XFEM to solve linear elastic fracture mechanics 
problems. The conventional expansion of the displacement field using a poly- 
nomial basis fails to capture the local behavior of the problem (e.g., steep 
stress gradients, material discontinuity). Hence, new functions are added 
to the conventional set of basis functions that contain information regard- 
ing the localized behavior. In general, the field variables are approximated 
by [351 [361 EH [381 [391 HO]: 

u'^(x) = 2, A^t(x)qi + enrichment functions (19) 

where A^i(x) are standard finite element shape functions, q/ are nodal vari- 
ables associated with node I. In the following, we briefly describe the stan- 
dard discretization of a plate using the field consistent Q4 plate element. 
The enriched field consistent Q4 element is then described. And finally, the 
discretized equations for the eigenvalue problem is given. In this section, 
only the essential details are given, interested readers are referred to recent 
review papers on XFEM ^T[ 1^ . A review on the implementation of the 
extended finite element is given in [13] and a detailed description of the state 
of the art on the simulation of cracks by partition of unity enriched methods 
is given in [H]. 
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3.1. Element Description 

The plate element employed here is a C'' continuous shear flexible field consis- 
tent element with five degrees of freedom («„, f o, Wq, O^, 6y) at four nodes in an 
4-noded quadrilateral (QUAD-4) element. If the interpolation functions for 
QUAD-4 are used directly to interpolate the five variables {uo,Vo,Wo,0x,9y) 
in deriving the shear strains and membrane strains, the element will lock 
and show oscillations in the shear and membrane stresses. The field consis- 
tency requires that the transverse shear strains and membrane strains must 
be interpolated in a consistent manner. Thus, the 6^ and 6y terms in the ex- 
pressions for shear strain Eg have to be consistent with the derivative of the 
field functions, Wo,x and Wo,y. This is achieved by using field redistributed 
substitute shape functions to interpolate those specific terms, which must 
be consistent as described in [211 ES]- This element is free from locking 
and has good convergence properties. For complete description of the el- 
ement, interested readers are referred to the literature [211 ES], where the 
element behavior is discussed in great detail. Since the element is based on 
the field consistency approach, exact integration is applied for calculating 
various strain energy terms. 

3.2. Enriched Q4 element 

Consider a mesh of field consistent Q4 elements and an independent crack 
geometry as shown in Figure ([3|). The following enriched approximation 
proposed by Dolbow et al., [45j for the plate displacements and the section 
rotations are used: 

{U\v\w'^) (X) = J2 <l>^{^){^^^,V^,Wi) + J^ U^)H{^){b] ,b] ,bj) + 

( 4 \ (20) 

^0.(x) 5^(c^„4„c^,)Gz(r,^) 

The section rotations in the shear terms are approximated by 

(^'' ^J) w = E <^^(^)(^-- ^..) + E 0.(x)i/(x)(6;^ 6,^'')+ 

/ 4 X (21) 

k&i^ \ 1=1 
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Figure 3: A typical FE mesh with an arbitrary crack. 'Squared' nodes are enriched 
with the heaviside function and 'circled' nodes with the near tip functions, which allows 
representing cracks independent of the background mesh. 



11 



In Equations (120]) and (1^ . (mj, Vi,Wi, 9^^, 9yJ are the nodal unknown vectors 
associated with the continuous part of the finite element solution, bj is the 
nodal enriched degree of freedom vector associated with the Heaviside (dis- 
continuous) function, and Cki is the nodal enriched degree of freedom vector 
associated with the elastic asymptotic crack-tip functions. The asymptotic 
functions, Gi and Fi in Equations f l2U]) and fl^ are given by ([15]): 

Gz(r,^) ^ l^sin (^0 , ^cos (^0 , ^sin (^^ , ^cos (j^ | , (22a) 

Fi{r, 0) = yf jsin f , cos (^j , sin (^j ^i^ (0) , cos (^j sin (6) | .(22b) 

Here (r, 6') are polar coordinates in the local coordinate system with the 
origin at the crack tip. The functions described here to recover the singular 
fields around the crack tip were originally proposed for isotropic plates |15] . 
As we are interested in the global behavior of the cracked FGM plate, we 
propose to use the same enrichment functions. The role of these enrichment 
functions is to aid in representing the discontinuous surface independent of 
the mesh. 

In Equations (!20!) and (I2T1) . 3\|^"=™ is the set of all nodes in the mesh; Jsf^ is 
the set of nodes whose shape function support is cut by the crack interior 
(squared nodes in Figure ([3])) and 3\f^ is the set of nodes whose shape function 
support is cut by the crack tip (circled nodes in Figure ([3])). For any node 
in y^, the support of the nodal shape function is fully cut into two disjoint 
pieces by the crack. If for a certain node tt,,, one of the two pieces is very 
small compared to the other, then the generalized Heaviside function used 
for the enrichment is almost a constant over the support, leading to an ill- 
conditioned matrix [46J. Therefore, in this case, the node rii is removed from 
the set 'N^. The area-criterion for the nodal inclusion in 'N^ is as follows: 
let the area above the crack is A"-^ and the area below the crack is A'^f and 
A^ = A"^ -f J^^. If either of the two ratios, A"^ j A^ or J^^ I A^ is below a 
prescribed tolerance, the node is removed from the set Ji'^ . A tolerance of 
10~^ is used in the computations. 

3.3. Discretized equations for enriched Q4 plate element 
Now, applying the displacement field approximated by Equations fl20|) - fl2Tl) 
in Equation flTll) and Equation flT6|) . one gets the modified Lagrange's equa- 
tions of motion in discretized form as: 
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(k - u^m) S = 0, 
where the element stiffness matrix is given by: 
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where B^^^ and 'Qenr are the standard and enriched part of the strain-displacement 
matrix, respectively and D is the material matrix. The element mass matrix 
is given by: 



M^ 



uu 



ua 



Ml, M^ 



KtdP'^std KtdP'^e 
^l^.p-Nstd KnrP'^e 



dQe- (25) 



where in deriving the above element mass matrix, the plate displacements 
and the section rotations given by Equations (12(1 - (1^ are used. In solving 
for the eigenvalues, the QR algorithm, based on the QR decomposition is 
used 



4. Numerical results 

In this section, we present the natural frequencies of a cracked function- 
ally graded material plates using the extended Q4 formulation. We consider 
both square and rectangular plates with simply supported, cantilevered and 
clamped boundary conditions. In all cases, we present the non dimensional- 
ized free flexural frequencies as, unless specified otherwise: 



Q 




(26) 



where E^., Vc are the Young's moulus and Poisson's ratio of the ceramic ma- 
terial, and pc is the mass density. In order to be consistent with the existing 
literature, properties of the ceramic are used for normalization. The effect of 
plate thickness a//i, aspect ratio 6/a, crack length d/a, crack orientation 9, 
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location of the crack, multiple cracks and boundary condition on the natu- 
ral frequencies are studied. Based on progressive mesh refinement, a 34 x 34 
structured mesh is found to be adequate to model the full plate for the present 
analysis. The material properties used for the FGM components are listed 
in Table [1] 

Before proceeding with the detailed study on the effect of different param- 
eters on the natural frequency, the formulation developed herein is vali- 
dated against available results pertaining to the linear frequencies of cracked 
isotropic and functionally graded material plates with different boundary 
conditions. The computed frequencies for cracked isotropic simply supported 
rectangular plate is given in Table [31 Tables S] and O gives a comparison of 
computed frequencies for simply supported square plate with a side crack 
and cantilevered plate with a side crack, respectively. It can be seen that 
the numerical results from the present formulation are found to be in good 
agreement with the existing solutions. 

The FGM plate considered here consists of silicon nitride (SI3N4) and stain- 
less steel (SUS304). The material is considered to be temperature dependent 
and the temperature coefficients corresponding to SI3N4/SUS304 are listed 
in Table |2] jlHl [31]. The mass density (p) and thermal conductivity (K) 
are: pc=2370 kg/m^, ^^=9.19 W/mK for SI3N4 and pm = 8166 kg/m^, K^ 
= 12.04 W/mK for SUS304. Poisson's ratio z/ is assumed to be constant 
and taken as 0.28 for the current study |1H1 [IH]- Here, the modified shear 
correction factor obtained based on energy equivalence principle as outlined 
in [50] is used. The boundary conditions for simply supported and clamped 
cases are (see Figure (jl])): 
Simply supported boundary condition: 



Uo = Wo = Oy = on X = 0, a 

Vo = Wo = Ox = ony = 0,b (27) 



Clamped boundary condition: 



Uo = Wo = Oy = Vo = Ox = on X = 0, a 



Uo = Wo = Oy = Vo = Ox = on y = 0,b (28) 
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Table 1: Material properties of the FGM components. ^Ref [2T]. * [311148] 



Material 




Properties 




E(GPa) 


V p 


(Kg/m3) 


Aluminum (Al)^ 


70.0 


0.30 


2702 


Alumina (Al203)^ 


380.0 


0.30 


3800 


Zirconia (Zr02)^ 


200.0 


0.30 


5700 


Steel (SUS304)* 


201.04 


0.28 


8166 


Silicon Nitride (Si3N4)* 


348.43 


0.28 


2370 



Table 2: Temperature dependent coefficient for material SI3N4/SUS3O4, Ref [511 H5] . 

Material Property P^ P_i Pi P2 P3 

E(Pa) 348.43e9 0.0 -3.070e-^ 2.160e-^ -8.9466"^^ 
SI N 

a (1/K) 5.8723e-6 0.0 9.095e-^ 0.0 0.0 

E(Pa) 201.04e9 0.0 3.079e-^ -6.534e-^ 0.0 
SUS304 

a (1/K) 12.330e-6 0.0 8.086e-4 0.0 0.0 



Table 3: Comparison of frequency parameters uj{h'^ /h)^ pd Ec for a simply supported 
homogeneous rectangular thin plate with a horizontal crack (a/6 — 2,b/h — 100, Cy/b = 
0.5, d/a = 0.5, 9 = 0). 

mode Ref [15] Ref M Ref ED Present 



1 


3.050 


3.053 


3.047 


3.055 


2 


5.507 


5.506 


5.503 


5.508 


3 


5.570 


5.570 


5.557 


5.665 


4 


9.336 


9.336 


9.329 


9.382 


5 


12.760 


12.780 


12.760 


12.861 



15 



Table 4: Non-dimensionalized natural frequency for a simply supported square AI/AI2O3 
plate with a side crack {a/b — 1, a/h — 10). Crack length d/a = 0.5. 



gradient 


Mode 1 


Mode 2 


Mode 3 


index, n 


Ref [21 


Present 


Ref [21 


Present 


Ref 21] 


Present 





5.379 


5.387 


11.450 


11.419 


13.320 


13.359 


0.2 


5.001 


5.028 


10.680 


10.659 


12.410 


12.437 


1 


4.122 


4.122 


8.856 


8.526 


10.250 


10.285 


5 


3.511 


3.626 


7.379 


7.415 


8.621 


8.566 


10 


3.388 


3.409 


7.062 


7.059 


8.289 


8.221 



Table 5: Fundamental frequency Lub"^ /h^J pd Ec for cantilevered square Al/Zr02 FGM 
plates with horizontal size crack [b/h = 10,Cj,/6 = 0.5, d/a = 0.5). 

a/h Mode gradient index, n 



0.2 1 5 10 



Ref [21J 1.0380 1.0080 0.9549 0.9743 0.9722 

Present 1.0380 1.0075 0.9546 0.9748 0.9722 

Ref |21j 1.7330 1.6840 1.5970 1.6210 1.6170 

Present 1.7329 1.6834 1.5964 1.6242 1.6194 

Ref [21] 4.8100 4.6790 4.4410 4.4760 4.4620 

Present 4.8231 4.6890 4.4410 4.4955 4.4845 
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4.I. Plate with a center crack 

Consider a plate of uniform thickness, h and with length and width as a and 
b, respectively. Figure (jl]) shows a plate with all edges simply supported with 
a center crack of length c. 



hy 



1 

1 i 


d 


1 


1 




Cy 




1 V 


1 



X 



Figure 4: Simply supported plate with a center crack. 

Effect of crack length, crack orientation and gradient index 

The influence of the crack length d/a, crack orientation 9 and gradient index 
n on the fundamental frequency for a simply supported square FGM plate 
with thickness a/h =10 is shown in Tables [6] and [3 It is observed that as the 
crack length increases, the frequency decreases. This is due to the fact that 
increasing the crack length increases local flexibility and thus decreases the 
frequency. Also, with increase in gradient index n, the frequency decreases. 
This is because of the stiffness degradation due to increase in metallic volume 
fraction. It can be seen that the combined effect of increasing the crack 
length and gradient index is to lower the fundamental frequency. Further, it 
is observed that the frequency is lowest for a crack orientation 6 = 45°. The 
frequency values tend to be symmetric with respect to a crack orientation 
6 = 45°. This is also shown in Figure ([5]) for gradient index n = 5 and crack 
length d/a = 0.8. 

Effect of crack location 

Next, the influence of the crack location on the natural frequency of a square 
plate with thickness, a/h = 10 and crack length, d/a = 0.2 is studied. The 
results are presented in Figure (Q. It is observed that the natural frequency 
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Table 6: Fundamental frequency uj}? jh^ p^jEc. 
square plate, ^denotes change in trend. 



for simply supported Si3N4/SUS304 FGM 



gradient Crack 



Crack length, dja. 



index, n 


orientation, d 





0.4 


0.6 


0.8 







5.5346 


5.0502 


4.7526 


4.5636 




10 


5.5346 


5.0453 


4.7386 


4.5337 




20 


5.5346 


5.0379 


4.7043 


4.4509 




30 


5.5346 


5.0278 


4.6640 


4.3528 




40 


5.5346 


5.0207 


4.6370 


4.2849 





45t 


5.5346 


5.0173 


4.6342 


4.2754 




50 


5.5346 


5.0204 


4.6370 


4.2849 




60 


5.5346 


5.0278 


4.6640 


4.3528 




70 


5.5346 


5.0380 


4.7043 


4.4509 




80 


5.5346 


5.0453 


4.7384 


4.5337 




90 


5.5346 


5.0503 


4.7527 


4.5636 







3.3376 


3.0452 


2.8657 


2.7518 




10 


3.3376 


3.0422 


2.8571 


2.7337 




20 


3.3376 


3.0376 


2.8362 


2.6833 




30 


3.3376 


3.0315 


2.8117 


2.6237 




40 


3.3376 


3.0271 


2.7953 


2.5825 


1 


45t 


3.3376 


3.0252 


2.7936 


2.5767 




50 


3.3376 


3.0270 


2.7953 


2.5824 




60 


3.3376 


3.0315 


2.8117 


2.6237 




70 


3.3376 


3.0377 


2.8363 


2.6833 




80 


3.3376 


3.0422 


2.8571 


2.7337 




90 


3.3376 


3.0452 


2.8657 


2.7518 



Table 7: Fundamental frequency uj}? jh^ p^jEc. 
square plate. ^ denotes change in trend. 



for simply supported Si3N4/SUS304 FGM 



gradient Crack 



Crack length, dja. 



index, n 


orientation, d 





0.4 


0.6 


0.8 







3.0016 


2.7383 


2.5769 


2.4747 




10 


3.0016 


2.7356 


2.5692 


2.4583 




20 


3.0016 


2.7315 


2.5504 


2.4130 




30 


3.0016 


2.7259 


2.5283 


2.3594 




40 


3.0016 


2.7220 


2.5136 


2.3223 


2 


45t 


3.0016 


2.7202 


2.5120 


2.3170 




50 


3.0016 


2.7219 


2.5135 


2.3222 




60 


3.0016 


2.7259 


2.5283 


2.3594 




70 


3.0016 


2.7315 


2.5504 


2.4130 




80 


3.0016 


2.7356 


2.5692 


2.4583 




90 


3.0016 


2.7383 


2.5770 


2.4747 







2.7221 


2.4833 


2.3371 


2.2445 




10 


2.7221 


2.4809 


2.3302 


2.2297 




20 


2.7221 


2.4772 


2.3131 


2.1887 




30 


2.7221 


2.4722 


2.2932 


2.1402 




40 


2.7221 


2.4686 


2.2798 


2.1067 


5 


45t 


2.7221 


2.4670 


2.2785 


2.1019 




50 


2.7221 


2.4685 


2.2798 


2.1066 




60 


2.7221 


2.4722 


2.2932 


2.1402 




70 


2.7221 


2.4772 


2.3132 


2.1887 




80 


2.7221 


2.4809 


2.3301 


2.2297 




90 


2.7221 


2.4833 


2.3371 


2.2445 
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10 20 30 40 50 60 70 80 90 

Crack orientation, 8 



Figure 5: Varition of fundamental frequency with orientation of the crack for a simply 
suppored square FGM plate, a/h = 10. 



of the plate monotonically decreases as the crack moves along the edges and 
towards the center of these edges. The natural frequency of the plate is 
maximum when the damage is situated at the corner. As the crack moves 
along the center lines of the plate from the edges and towards the center 
of the plate, the natural frequency increases up to a certain distance and 
then decreases. When the crack is situated at the center of the plate, the 
frequency is minimum. 

Plate with multiple cracks 

Figure ([7]) shows a plate with two cracks with lengths ai and 02 and with ori- 
entations 9i and 62 they subtend with the horizontal. The effect of respective 
crack orientation on the fundamental frequency for a simply supported FGM 
plate is numerically studied. The horizontal and vertical separation (see 
Figure ([7])) between the crack tips is set to a constant value, H =0.2 and 
V =0.1, respectively. Table |8] shows the variation of fundamental frequency 
for plate with gradient index n =5 and crack lengths Oi = 02 = 0.2 as a 
function of crack orientations. Figure ([H]) shows the variation of frequency as 
a function of orientation of one crack for different orientations of the second 
crack. It can be seen that with increase in crack orientation, the frequency 
initially decreases until it reaches a mimimum at 6' = 45°. With further 
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2.69 



2.695 



2.7 



2.705 



2.71 



2.715 




Along x-direction 0.2 



Along y-directlon 



(a) 




-0.1 0.1 

Along x-dlrectlon 



(b) 

Figure 6: Variation of fundamental frequency as a function of crack position for a simply 
supported square FGM plate. The crack orientation, 9 is taken to be 0, i.e., horizontal 
crack. 
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increase in crack orientation, the frequency increases and reaches maximum 
at 9i = 62 =90°. The frequency value at ^ = 0° is not the same as that at 
6 = 90° as we would expect. This is because when 9i = 62 = 90°, the crack is 
located away from the center of the plate and crack disturbs the mode shape 
shghtly. 




Figure 7: Plate with multiple cracks: geometry 
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Table 8: Fundamental frequencj 
denotes change in trend. 


' c^byh^. 


~/Ec for simply supported Si3N4/SUS304 FGM square plate, gradient index, n ■ 


=5. 




Crack 










Crack orientation, 62- 










Orientation, 6i 





10 


20 


30 


40 


45 


50t 


60 


70 


80 


90 





2.6149 


2.6122 


2.6077 


2.6034 


2.5997 


2.5988 


2.5981 


2.5993 


2.6023 


2.6066 


2.6098 


10 


2.6122 


2.6094 


2.6049 


2.6006 


2.5969 


2.5960 


2.5953 


2.5964 


2.5994 


2.6037 


2.6068 


20 


2.6077 


2.6049 


2.6004 


2.5961 


2.5924 


2.5914 


2.5907 


2.5919 


2.5948 


2.5991 


2.6022 


to 30 

00 


2.6034 


2.6006 


2.5961 


2.5918 


2.5881 


2.5871 


2.5864 


2.5876 


2.5905 


2.5948 


2.5979 


40 


2.5997 


2.5969 


2.5924 


2.5881 


2.5844 


2.5835 


2.5827 


2.5839 


2.5868 


2.5910 


2.5941 


45 


2.5987 


2.5959 


2.5914 


2.5871 


2.5834 


2.5825 


2.5818 


2.5829 


2.5858 


2.5901 


2.5931 


50t 


2.5980 


2.5952 


2.5907 


2.5864 


2.5827 


2.5818 


2.5811 


2.5822 


2.5851 


2.5893 


2.5924 


60 


2.5992 


2.5964 


2.5918 


2.5876 


2.5839 


2.5829 


2.5822 


2.5834 


2.5863 


2.5905 


2.5935 


70 


2.6022 


2.5993 


2.5948 


2.5905 


2.5868 


2.5858 


2.5851 


2.5863 


2.5891 


2.5933 


2.5963 


80 


2.6065 


2.6036 


2.5990 


2.5947 


2.5910 


2.5900 


2.5893 


2.5904 


2.5933 


2.5974 


2.6004 


90 


2.6098 


2.6069 


2.6022 


2.5979 


2.5942 


2.5932 


2.5924 


2.5935 


2.5964 


2.6005 


2.6035 



Effect of aspect ratio, thickness and boundary conditions 

The influence of plate aspect ratio b/a, plate thickness a/h and boundary 
condition on a cracked FGM plate with a horizontal center crack is shown in 
Table M Two types of boundary conditions are studied: Simply supported 
(SS) and Clamped condition (CC). For a given crack length and for a given 
crack location, decreasing the plate thickness and increasing the plate aspect 
ratio, increases the frequency. The increase in stiffness is the cause for in- 
crease in frequency when the boundary condition is changed from SS to CC 
for a fixed aspect ratio and plate thickness. 



Table 9: Effect of plate aspect ratio b/a, plate thickness a/h and boundary condition 
on fundamental frequency wb"^ /hyj pc/Ec for Si3N4/SUS304 FGM plate with horizontal 
center crack {cy/b — 0.5, d/a — 0.5). 'Simply Supported, *Clamped Support. 



b/a 


a/h 


Mode 1 
SSt CC* 


Mode 2 




sst 


CC* 




10 


1.1205 


2.2202 


2.1586 


2.8588 


0.5 


20 


1.1974 


2.5043 


2.5482 


3.4621 




100 


1.2625 


2.6748 


2.6593 


4.0903 




10 


2.4051 


4.1624 


5.2792 


6.6286 


1 


20 


2.4831 


4.4592 


5.8338 


7.6464 




100 


2.5473 


4.6311 


6.2765 


8.5774 




10 


6.6864 


12.7513 


10.5295 


15.9115 


2 


20 


6.8101 


13.4551 


10.8647 


17.0301 




100 


6.8847 


13.7540 


11.0392 


17.6030 



4-2. Plate with side crack 

Consider a plate of uniform thickness, h, with length and width as a and 
b, respectively. Figure ([9]) shows a cantilevered plate with a side crack of 
length d located at a distance of Cy from the x-axis and at an angle 6 with 
respect to the x— axis. The influence of plate thickness, crack orientation 
and gradient index on the fundamental frequency is shown in Table [10] and 
in Figure ( ITOj) . With increase in gradient index n, the frequency decreases 
for all crack orientations and for different plate thickness. With increase in 
crack orientation from 6 = —60° to ^ = 60°, the frequency initially decreases 
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10 20 30 40 50 

Crack orientation, i 



(a) Effect of crack orientation, 



80 


-.^ 


70 


' T.?' 


60 




50 


T- 


40 


'^■e 


30 




20 




10 


^^•V 




10 20 30 40 50 60 70 80 90 



(b) Effect of gradient index, n 



Figure 8: Variation of fundamental frequency for a simply supported square Si3N4/SUS304 
FGM plate with a center crack as a function crack orientation with H =0.2, V =0.1 and 
a/h =10 (see Figure ©). 
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until 9 = —40° and then reaches the maximum when the crack is horizontaL 
And with further increase in the crack orientation, the frequency decreases 
and the plate's response is symmetric. This is because, when the crack 
is horizontal {6 = 0), the crack is aligned to the first mode shape and the 
response of the plate is similar to the cantilevered plate without a crack. The 
frequency of the plate without a crack is greater than a plate with a crack. 
Figures fITT]) and ( 1121) shows first two mode shapes for a cantilevered plate 
with and without a horizontal crack. As explained earlier the frequency and 
the first mode shape for a plate with and without a crack are very similar. 
For any other crack orientation, the mode shape would be influenced by the 
presence of the crack. Figures (1131) and fll4p shows the first mode shape of a 
cantilevered plate with a side crack with orientations 6 = ±40° and 6 = ±60°, 
respectively. 




Figure 9: Cantilevered plate with a side crack: geometry 



5. Conclusion 

Natural frequencies of cracked functionally graded material plate is studied 
using the extended finite element method. The formulation is based on first 
order shear deformation theory for plates and a four-noded field consistent 
enriched element is used. The material is assumed to be temperature de- 
pendent and graded in the thickness direction. Numerical experiments have 
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Table 10: Fundamental frequency iob"^ / h^/pjEl for cantilevered square plate 
Si3N4/SUS304 FGM plate with a side crack {cy/h = 0.5, d/a = 0.5) as a function of 
crack angle and gradient index, ^denotes change in trend. 



a/h crack 



gradient index, n 



angle, 9 





1 


2 


5 


10 


-60 


0.9859 


0.5918 


0.5322 


0.4838 


0.4610 


-50 


0.9840 


0.5906 


0.5312 


0.4829 


0.4601 


-40^ 


0.9838 


0.5905 


0.5311 


0.4828 


0.4600 


-30 


0.9862 


0.5919 


0.5323 


0.4840 


0.4611 


-20 


0.9900 


0.5943 


0.5344 


0.4859 


0.4630 


-10 


0.9936 


0.5964 


0.5364 


0.4876 


0.4646 


10 ot 


0.9951 


0.5973 


0.5372 


0.4884 


0.4653 


10 


0.9936 


0.5964 


0.5364 


0.4876 


0.4646 


20 


0.9900 


0.5943 


0.5344 


0.4859 


0.4630 


30 


0.9862 


0.5919 


0.5323 


0.4840 


0.4611 


40t 


0.9838 


0.5905 


0.5311 


0.4828 


0.4600 


50 


0.9840 


0.5906 


0.5312 


0.4829 


0.4601 


60 


0.9859 


0.5918 


0.5322 


0.4838 


0.4610 


-60 


0.9949 


0.5972 


0.5371 


0.4883 


0.4653 


-50 


0.9927 


0.5959 


0.5359 


0.4872 


0.4643 


-40^ 


0.9924 


0.5957 


0.5357 


0.4871 


0.4641 


-30 


0.9944 


0.5969 


0.5368 


0.4881 


0.4651 


-20 


0.9979 


0.5989 


0.5387 


0.4898 


0.4667 


-10 


1.0011 


0.6009 


0.5404 


0.4913 


0.4682 


20 Ot 


1.0024 


0.6016 


0.5411 


0.4920 


0.4688 


10 


1.0011 


0.6009 


0.5404 


0.4913 


0.4682 


20 


0.9979 


0.5989 


0.5387 


0.4898 


0.4667 


30 


0.9944 


0.5969 


0.5368 


0.4881 


0.4651 


40t 


0.9924 


0.5957 


0.5357 


0.4871 


0.4641 


50 


0.9927 


0.5959 


0.5359 


0.4872 


0.4643 


60 


0.9949 


0.5972 


0.5371 


0.4883 


0.4653 
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(a) Effect of aspect ratio, a/h 
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>0.56 



£ 0.54 



0.52 



0.5 



0.48^ 
-60 



— «— 


n=1 


-e - 


n=2 


-B- 


n=5 



(»- - £)_ _0- - © 



_ -Or 



. e - -G- -©. _ , 



-o- - e - -o — <> 



-40 



-20 20 

Cracl< orientation, 6 



40 



60 



(b) Effect of gradient index, n 

Figure fO: Variation of fundamental frequency for a cantilevered square Si3N4/SUS304 
FGM plate with a side crack as a function of crack orientation, d/a =0.5, Cj,/a=0.5. 
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(0. = 0.4885 




(a) Without crack, wi =0.4885 



(0. = 0.4883 




(b) With crack, cji =0.4883 

Figure 11: First Mode shape for a cantilevered plate with a side crack with 9 ~ Q,n 
5, Cy/b = 0.5, d/a = 0.5, a/h = 10, b/a = 1. 
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CO =1.1608 




(a) Without crack, W2 =1.1608 



i»„ = 0.8223 




(b) With crack, cjs =0.8223 



Figure 12: Second Mode shape for a cantilevered plate with a side crack with 9 
5, Cy/b = 0.5, d/a = 0.5, a/h = 10, b/a = 1. 



0,n 



30 



CO, = 0.4828 




(a) e = 40° 



(0. = 0.4828 




(b) e = -40° 

Figure 13: First Mode shape (wi — 0.4828) for a cantilevered plate with a side crack with 
n = 5, Cy/b = 0.5, dja = 0.5, ajh = 10, bja = 1. 
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(0. = 0.4838 




(a) e = 60° 



(0. = 0.4838 




(b) e ^ -60° 

Figure 14: First Mode shape (wi — 0.4838) for a cantilevered plate with a side crack with 
n = 5, Cyjb = 0.5, dja = 0.5, a//i = 10, bja = 1. 



32 



been conducted to bring out the effect of gradient index, crack length, crack 
orientation, crack location, boundary condition, plate aspect ratio and plate 
thickness on the natural frequency of the FGM plate. Also the influence 
of multiple cracks and their relative orientation on the natural frequency is 
studied. From the detailed numerical study, the following can be concluded: 

• Increasing the crack length decreases the natural frequency. This is due 
to the reduction in stiffness of the material structure. The frequency is 
lowest when the crack is located at the center of the plate. 

• Increasing gradient index n decreases the natural frequency. This is 
due to the increase in metallic volume fraction. 

• Decreasing the plate thickness a/h and increasing the aspect ratio h/a 
increases the frequency. 

• For a cantilevered plate with a side crack, horizontal crack has the 
maximum frequency and the trend changes at ^ = ±40°. 

• Crack orientation 9 = 45° has been observed to be a critical angle. 
At this crack orientation, the frequency changes its trend for a square 
plate. 

• Increasing the number of cracks, decreases the overall stiffness of the 
plate and thus decreasing the frequency. The frequency is lowest when 
both the cracks are oriented at ^ = 50°. 
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